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^— V \ Abstract. A brief introduction to theories of the gravitational field with a Lagrangian that is a 

fvq ■ function of the scalar curvature is given. The emphasis will be placed in formal developments, 

, ' while comparison to observation will be discussed in the chapter by S. Joras in this volume. 
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'^ INTRODUCTION 

Q^' The predictions of General Relativity (GR) are confirmed to an impressive degree by 

O^' observations in number of situations [1]. In spite of this fact, theories that differ from 

^^! GR either in the limit of low or high curvature have been intensively studied lately, 

and have a long tradition, starting with a paper by Weyl in 1918 [2]. Although Weyl's 

motivation was related to the unification of GR and Electrodynamics, the current revival 

of these theories is twofold. In the case of low curvature, the aim is to describe the 

C<^ \ accelerated expansion of the universe that follows from several observations [jsl] (when 

QQ ■ interpreted in the standard cosmological model [4]) [l|. Regarding the high-curvature 

regime, it is important to note that there is no observational evidence of the behaviour 

of the gravitational field for very large values of the curvature. This makes objects such 

(3 ' as black holes and neutron stars the ideal places to look for deviations from General 

Relativity in the strong regime. In fact, the Kerr solution is not unique in f{R) theories 

[JMJ. Consequently, any deviation from Kerr's spacetime in compact objects will be 

unequivocally signaling the need for changes in our description of strong gravity. The 

^ ■ task of understanding what kind of deviations can be expected, and their relation to 

observable quantities is of relevance in view of several developments that offer the 

prospect of observing properties of black holes in the vicinity of the horizon [12f|. 

In this short review we will be concerned with gravitational theories described by the 
action 

S= fd^Xy^fiR), (1) 



(N 
> 
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where g is the determinant of the metric ^^v and / is an arbitrary function of the 
curvature scalar R q The function / must satisfy certain constraints, some of which 



^ The possibility of describing the current accelerated expansion of the universe using f{R) theories was 
first discussed in Js]]. 

This choice is favoured by a theorem by Ostrogradski lll3ll over Lagrangians built with invariants 
obtained from the the Ricci and Riemann tensors. 



are necessary for the theory to be well-defined ab initio, and others to account for 
observational facts. Those in the first class will be discussed in this review, while those 
in the second class are presented in the chapter by S. Joras in this volume B We shall 
begin by reviewing in the next section some general features of this type of theory. 

THE THREE VERSIONS OF f{R) THEORIES 

We shall see in this section that f{R) theories can be classified in three different types, 
according to the role attibuted to the connection. In all of the versions, the equation for 
the energy-momentum conservation is valid, since the total (gravitational plus matter) 
action is diffeomorphism-invariant and gravity and matter are minimally coupled by 



hypothesis (see for instance [llll] . MM)- 



Metric version 

In the metric version of f{R) theories, the action 

S = ^ I d'^xV^fiR) +SM{gf,y,¥), (2) 

is varied with respect to g^v Here, 5m is the matter action, which is independent of the 
connection. The resultant equations of motion are of fourth order in the derivatives of 
the metric tensor: 

-^^^/iv - :^f{R)gfiv - [V^ Vv - g^vO\ — ^ = fcr^v, (3) 

where T^y is the energy-momentum of he matter fields, defined by 

2 55m 
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and the covariant derivative is defined using the usual Levi-Civita connection. Taking 
the trace, we obtain 

which is to be compared to R = —kT, the result in GR. 



^ There are several reviews that deal with different aspects of f{R) theories, see |Ia],|l3],|la],|l21, and also 
the recent book JTui . 



Equivalence with Brans-Dicke theory 

As shown for instance in [14], the gravitational part of the action given in Eqn.© is 
equivalent to the following action: 



d Xy 

where 



^-(*) 



(4) 






^ = f,x{x)' ^iid X = R' corresponding to a Brans-Dicke theory with (O — [|l6ll . Note 
that the absence of a kinetic term for the scalar field does not mean that it is non- 
dynamical: its evolution, due to the non-minimal coupling with the gravitational field, is 
given by the variation of the action wrt ^ : 

3n(^+2U{(j>)-^^ = KT. (6) 

Through a conformal transformation of the metric and a redefinition of the scalar field, 
the action given in Eqn.dH) can be taken to that of a scalar field minimally coupled 
with gravity, and with nonzero kinetic term and potential. These representations of f{R) 
theories show that there is a massive scalar degree of freedom, which manifests as a 



longitudinal mode in gravitational radiation (see for instance [|24|1 for the cosmological 
case). 

We close this section by stating that the equivalent representations are convenient 
since the associated equations of motion are of order two, but a word of caution is needed 
because sometimes the potential in Eqn.© is typically multivalued (see for instance 

I). It may be better to work directly in the original representation, as for instance 

1511 in the case of compact stars. 



[inDiis 



in Ref. 



Palatini version 

In this second type of f{R) theories, the metric and the connection are taken as 
independent fields, and the matter action 5m is independent of the connection. So the 
starting point is the action 

where ^ = g^^^M^v, and the corresponding Riemann tensor is constructed with a 
connection F a priori independent of the metric. 



From the variation of the action wrt the metric and F we get|j (see for instance t^) 

f'{^)^{nv) - ^^fmsnv = kT^v, (8) 

^xiv^fmsn=o, (9) 

where the prime denotes derivative wrt ^, and the barred covariant derivative is built 
with the connection F. GR is recovered by setting f{^) = ^ in these equations. Taking 
the trace of Eqn.® we obtain 

f'i^)^-2f{^) = KT, (10) 

which shows that in this case the relation between ^ and T is algebraic, hence no scalar 
mode is present. 

From Eqn.® , it follows that lH] 

rjv = fi^s^'" [Mf'i^)syo) + dy{f'mg,^) - dM'msnv)] ■ (11) 

Since this expression relates T to M and the metric, and ^ and T are in principle 
interchangeable through Eqn.(fTOl). the connection can be in principle expressed in terms 
of the matter fields and the metric. In other words, it is an auxiliary field. Li fact, Eqn.dH) 
can be rewritten as 

G^v = jjT^v-lg^vl^^-j]+j:j{y^yv-gnvD)f (12) 
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where the Einstein tensor and the covariant derivatives are built with the Levi-Civita 
connection, and M is expressed in terms of T using Eqn. (fTOl) . It follows that this version 
of f{R) theories can be interpreted as GR with a modified source. Perhaps the most 
important modification is that third order derivatives of the matter fields appear on the 
rhs of Eqn.(fT3l). As reported in [25], this feature may cause serious problems in static 
spherically symmetric solutions with a polytropic fluid with index 3/2 < 7 < 2 as a 
source. Note however that this result was challenged in the review [|26ll q 



In the case of GR, this method furnishes the same result as the metric case, but this is not the case in 
more general theories as discussed for example in {M\. 

^ It has also been claimed that the Cauchy problem is not well-posed for the Palatini version of f{R) 
theories j^. For an updated discussion see ||27|] . 



Metric-affine theories 

In this case, the matter action depends of the connection, which is a priori independent 
of the metric. The action is given by 

S=^j d^x^^gf{M)+SM{giiv,Tlv.W)- (13) 

Depending on the matter fields, the theory may display non-propagating torsion and 
non-metricity (see [|90 and [ | 36l] for details). 

NONMINIMAL COUPLING 

Metric f{R) theories have been generalized by allowing a nonminimal coupling between 
the curvature and the matter Lagrangian, with action given by 

S = J[\MR) + [1 + A/2(^)]i^n} V^d\ (14) 

where /i and fj are arbitrary functions of R, and A is a constant. A particular case of 
this action was considered in [28] in the context of the accelerated expansion of the 
universe. Later, it was shown in [29] u that this type of theory leads to a modification of 
the conservation law of the matter energy-momentum tensor, namely 

V^# = T-^/^ Lv^m - #1 V^R. (15) 



1+A/: 
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The presence of a nonzero rhs leads to non-geodesic motion, and it was suggested in 



U% that this may be related to MOND. 
A more general type of theories was proposed in [|30D, with action given by 



S = lf{R,^m)v^d\ (16) 

where / is an arbitrary function of R and of the matter Lagrangian. As in the previous 
case, an extra force, perpendicular to the 4-velocity, accelerates the particles. To close 
this section, let me mention that all the versions of f{R) theories mentioned above can 
be described in a unified framework in the so-called C-theory [|42|] . 

ASSORTED APPLICATIONS 

f{R) theories have been used to describe different aspects of relativistic astrophysics and 
cosmology. Since the low curvature limit, which has been studied primarily to explain 
the accelerated expansion of the universe, is discussed in the chapter by S. Joras in 
this volume, only one example will be given here in this regime. Afterwards, some 
applications in the strong-curvature regime will be discussed. 



^ See also H 



Low curvature 

In the case of the k — Friedmann-Lemaitre-Robertson- Walker metric, the EOM ([3]) 
can be written as 

p=-/x-{+3r^, (17) 

2 a 



P = -^iRn+R) + {-f"(R-^yrR'. 



(18) 
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Let us remark that it is safe to asssume that most of the current matter content of the 
universe (assumed here to be normal matter, as opposed to dark energy) is pressureless. 
This matter must satisfy the conditions po > and po = 0. where the subindex means 
that the quantity is evaluated today. Using Eqns.(fT7l) and (fTSi) . we shall rewrite these 
conditions in terms of following kinematical parameters: the Hubble and deceleration 
parameters, the jerk, and the snap, respectively given by [ 431 

a 
a 

While the current value of the first two parameters is relatively well-established today, 
the value of jo is not known with high precision, and no acceptable value of ^o has been 
reported yet [441 . By writing po > in terms of the kinematical parameters we get 

3qoHU-^-im^f;^Uo-qo-2)>0. (19) 

This inequality gives a relation that the parameters and the derivatives of a given f{R) 
must satisfy today and, as shown in [451, it limits the possible values of the parameters 
of a given theory. Notice that Eqn.([T8l) involves the snap (through R). If we had a 
measurement of sq, we could use the equation po = to obtain another constraint on 
f{R). Since this is not the case, we shall express po = in such a way that it gives a 
forecast fo the possible current values of the snap for a given f{R) : 

rf rill r 

'0 = T^iqo-2)+6Hi^-^i-qo + Jo-2f-[qoiqo + 6)+2{l+Jo)]~j^^ 
0/io/o Jo ^^^ h 

(20) 



Strong curvature 

The possibilit y o f nonsingular cosmological solutions in f{R) theories has been 



considered in [[31[1 and [321. In the latter article, a necessary condition for a bounce 
to occur in a Friedmann-Lematre-Robertson-Walker setting was obtained, and it is 
given by 

^ _ _Po ^ (-91^ 

«0 Jo 2JQ 



with 

i?o = 6f^ + 4Y (22) 

and the subindex b means that the quantity is evaluated at the bounce. Contrary to 
the case of GR, a bounce may occur for any value of K. 
It was shown in [|33ll that the theory given by 

has an inflationary solution given by 

H^H,-—{t~t,), 

D 



a ~ a,- exp 



■2 



M 

Hi{t-ti)-—{t-t, 
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where tj marks the beginning of the inflationary epoch. Several features of this 
model have been studied in detail (see references in [8]). The results of WMAP 
constraint M ^ lO^'^ GeV, and the spectral index for this model is n<yg ^ 0.964, 
which is in the range allowed by WMAP 5-year constraint. The tensor to scalar 
ratio r also satisfies the current observational bound, but is different from that of 
the chaotic inflation model. Hence, future observations such as the Planck satellite 
may be able to discriminate between these two models. 

Compact stars have been repeatedly studied for a number of f{R) theories, either 
in the conformal representation [ | 37i1 . or directly in the fourth-order version (see for 
instance [15]). 

Regarding black holes in f{R) theories, it was shown in |34l ] that an observational 
verification of the Kerr solution for an astrophysical object cannot be used in 
distinguishing between GR and f{R) theories. Hence, the observation of deviations 
from the Kerr spacetime may point to changes in our understanding of gravitation. 



Other features of black holes in f{R) theories have been analyzed in [|35|1 



CONCLUSIONS 

In this short review, I intended to show that several aspects of f(R) theories (in its various 
representations) have been extensively studied in the literature. There are many other 



aspects that I had to leave aside such as "good propagation" (i.e. absence of shocks) [|38|] . 
the loop representation ll39ll . and the Hamiltonian representation [40]. Although the low- 
curvature regime and its consequences has attracted a lot of attention due to its possible 
relevance in Cosmology, the high-curvature regime is also of interest independently of 
the low-curvature regime, and the consequences of a modification in such a regime 



may be testable in the near future (as in the inflationary model in [|33|] . and through the 



observation of electromagnetic [12] and gravitational waves [|41ll in the case of compact 
objects). 
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